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Abstract. In this study we explore the magnetic mechanism of hypernovae and relativistic jets of long duration gamma ray 
bursts within the collapsar scenario. This is an extension of our earlier work[T|. We track the collapse of massive rotating 
stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a 
realistic equation of state and takes into account the cooling associated with emission of neutrinos and the energy losses 
due to dissociation of nuclei. The neutrino heating is not included. We describe solutions with different black hole rotation, 
mass accretion rate, and strength of progenitor's magnetic field. Some of them exhibits strong explosions driven by Poynting- 
dominated jets with power up to 12 x 10 51 erg s _1 . These jets originate from the black hole and powered via the Blandford- 
Znajek mechanism. A provisional criterion for explosion is derived. A number of simulation movies can be downloaded from 
http : //www. maths . leeds . ac .uk~ serguei/research/movies/anim.html 
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INTRODUCTION 

The most popular model of central engines of long- 
soft Gamma Ray Burst (GRB) is based on the "failed 
supernova" scenario of stellar collapse, or "collapsar", 
where the core of a rapidly rotating progenitor star forms 
a black hole (BH) surrounded by an accretion disk 0. 
A similar configuration can be produced via inspiraling 
of a BH or a neutron star into the companion star during 
the common envelope phase of close binary (3j (4) . The 
gravitational energy released in the disk can be very 
large, more than sufficient to stop the collapse of outer 
layers and drive GRB outflows, presumably in the polar 
direction where density is much lower 0. GRB jets 
can be powered via heating due to annihilation of disk 
neutrinos |5 ] and via magnetic braking (SHE). 

The high angular momentum assumed for the stellar 
core in the collapsar model also implies high rotation 
rates of the central black hole and the possibility of pow- 
ering the GRB jets via the Blandford-Znajek mechanism 
l9l . The rotational energy of black hole with mass Mh = 
2Mq and spin a = 0.9 is enormous, E rot ^ 5 x 10 53 erg, 
which is fifty times higher than the rotational energy of a 
millisecond pulsar [10 ] and well above what is needed to 
explain GRBs and associated hypernovae. Applying the 
Blandford-Znajek formula to the collapsar problem we 
obtain the following estimate of power 

E BZ = 3.6 x lO^M^M^eigs" 1 , (!) 

where M 2 = M h /2M Q , ¥ 27 = ^/lO^Gcm 2 is the mag- 
netic flux penetrating BH, f(a) = a 2 /(I + \f\ —a 1 ) 1 , 



and it is assumed that BH magnetosphere rotates with 
angular velocity £1 = 0.5^, thus delivering maximum 
power. These estimates show that the magnetic braking 
of BH could fully account for the energetics of GRBs. 
However, the result £2 = 0.5^ is obtained for force-free 
magnetospheres where one can ignore the inertia of mat- 
ter and it is not clear when and how such conditions can 
develop in the collapsar model. For high mass-loading 
in the magnetosphere one could expect the magneto- 
spheric rotation to be dictated by the angular momentum 
of matter rather then by the rotation of BH itself. Careful 
GRMHD simulations are needed in order to resolve this 
issue. Previous studies, with simplified setting and micro 
physics, have already demonstrated that this problem can 
be addressed with modern computational tools (e.g. (HI) 
and invite to explore progressively more realistic models. 
Here we present numerical solutions for models with re- 
alistic EOS and effects of neutrino cooling. The neutrino 
transport, however, is rather complicated, particularly in 
the curved space of BH, and for this reason we have not 
included the heating due to neutrino annihilation. In re- 
ality, both the neutrino and magnetic mechanisms could 
operate hand in hand and we plan to include neutrino 
heating in future studies. The results for one particular 
numerical model have been presented in Q); here we out- 
line a wider study. The effects of neutrino cooling on disc 
accretion have been described in lfT2l . 



SIMULATION SETUP 

Since our code can deal only with time-independent 
space-time we are forced to start from the point where the 
central BH of mass Mh = 3M has already been formed 
inside the collapsing star. The rotation rate of BH is 
rather uncertain and depends on the angular momentum 
distribution in the progenitor. Here we use either the 
rather optimistic value of a = 0.9 or the very pessimistic 
value of a = 0.0. The latter was introduced to fully 
eliminate the Blandford-Znajek effect and to see if the 
magnetic braking of accretion disk alone could make the 
star to explode. 

The main details of our numerical method and various 
test simulations are described in fT3l . The only really 
new feature here is the introduction of HLL-solver which 
is activated when our linear Riemann solver fails, this 
occurs in regions of relativistically high magnetization. 
We have also found that total switch to the HLL-solver 
significantly degrades numerical solutions via increasing 
numerical diffusion and dissipation (see also fT4ll ). 

The gravitational attraction of BH is introduced via 
Kerr metric in Kerr-Schild coordinates, {£,0,r, 0}. The 
two-dimensional computational domain is (ro < r < 
n) x (0 < 6 < it), where r = (1 + 0.5\/l -a 2 ) r g and 
ri 5700 r g = 25000 km. 

Notice that the inner boundary is inside of the outer 
event horizon - we can do this because the horizon coor- 
dinate singularity is eradicated in the Kerr-Schild coor- 
dinates. The total mass within the domain is small com- 
pared to the mass of BH ( only about 7%) that allows us 
to ignore its self-gravity. The grid is uniform in 6 where 
it has 180 cells and almost uniform in log(r) where it 
has 450 cells, the linear cell size being the same in both 
directions. 

A realistic equation of state (EOS) is introduced using 
the EOS code HELlvQ [ 1511 . The neutrino cooling is 
computed using the interpolation formulas given in (TSJ 
ITTl H8l . The photo-disintegration of nuclei is included 
via modifying EOS as in lfT9l . The equation for mass 
fraction of free nucleons is adopted from [20]. 

The collapsing star is described by the free-fall model 
of EH . The collapse time for a core of radius r c = 10 9 cm 
and mass 3.OM is t c = Vlr c /3vff ~ 0.5s. Since, in 
this study we explore the possibility of explosion soon 
after the core collapse, we start simulations 1 second 
after the onset on collapse. Because GRBs are currently 
associated with very massive progenitors we set the mass 
parameter C\ of Bethe's model (see Q ) to either C\ = 
3 or Ci =9. These correspond to the accretion rates 



the code can be downloaded from the website 
i://www.cococubed.com/code_pages/eos.shtml 
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FIGURE 1. Model C\ = 9, B = 3 x 10 10 G at time t = 
0.24s, soon before the explosion. Colour shows the ratio of 
gas and magnetic pressure, log 10 (P/P m ), contours the magnetic 
field lines and arrows the velocity field. Notice the outflow 
beginning to develop in the highly magnetized "blue stripe" 
stretching at the angle of 45° from the black hole. 

of M ~ O.166M s _1 and M ~ O.5M s _1 respectively. 
The angular momentum distribution describes a solid 
body rotation within the cylindrical 77 = 6300km and is 
constant with the value of /o = 10 17 cm 2 s -1 further out, 
similar to that used in [5]. These initial conditions are 
also used to fix the flow variables in the ghost cells of the 
outer boundary. 

The magnetic field distribution is that of a uniformly 
magnetized sphere in vacuum, the radius of this sphere 
r\ = 4500km. Inside the sphere the magnetic field 
strength is B = 3 x 10 10 G, 10 10 G or 3 x 10 9 G. 

COMPUTER SIMULATIONS 

At the beginning of simulations the angular momentum 
of gas near BH is less than that of the last stable orbit, 
l\ so , and as a result it falls straight into the BH. The most 
rapidly varying parameter at this stage is the magnetic 
flux threading the BH which grows linearly with time. 
Later, when matter with / ~ l\ so reaches the BH it forms 
a thin inviscus accretion disk in the equatorial plane |22|. 
Finally, when the angular momentum of gas approaching 
the BH in the equatorial direction well exceeds l\ so the 
centrifugal force halts its infall and a turbulent accretion 
disk develops around the BH. At the same time a strong 
accretion shock lifts from its surface. After passing the 




accretion shock the low angular momentum plasma of 
polar regions keeps falling straight into the BH whereas 
the high angular momentum plasma coming at interme- 
diate directions fills the low density volume above and 
below the disk. Strong differential rotation within this 
bubble results in strong amplification of azimuthal mag- 
netic field. Further evolution depends of the amount of 
magnetic flux accumulated by the BH. 

In the model with C\ = 9 and B = 3 x 10 10 G we ob- 
served strong explosion very soon after formation of the 
accretion disk. The BH is a key factor in this explosion 
pumping the electromagnetic energy into the bubble at 
the rate of ~ 12 x 10 51 ergs -1 . Prior to the explosion we 
see mass unloading of some magnetic field lines connect- 
ing the bubble with BH (see fig[T]). This seems to be pro- 
moted by the fact plasma entering the accretion shock 
along these field lines is diverted towards the equatorial 
plane instead of entering the bubble and falling into the 
BH. Similar evolution was shown by the model described 
in 

In the model with C\ = 3 and Bq = 10 10 G we observe 
several pulsations of the accretion shock, gradual expan- 
sions followed by relatively rapid contractions, prior to 
the explosion which is significantly less energetic. How- 
ever, once the explosion sets up it shows the same generic 
features (see fig(2]). We observe two well defined polar 
jets surrounded by magnetic cocoons of high pressure 
and low density. The magnetic pressure of these cocoons, 
which have been inflated by the jets, exceeds by more 
than six orders of magnitude the magnetic pressure in 
the collapsing star. These over-pressured cocoons drive a 
blast wave into the star. The mean propagation speed of 
the blast in the polar direction v s ~ 0. lc. In the vicinity of 
BH the solution shows the same key features as found in 
the previous studies of thick disks around BHs - main 



disk, its dynamic corona, and magnetically-dominated 
funnel inside which two Pointing-dominated jets are pro- 
duced bv the BH t23ll24l[25llTl[26l[T2ll. Integrating over 
the volume of the blast wave we can measure the rate of 
energy supply in the explosion, E ~ 0.4 x 10 51 ergs -1 . 
The direct measurement of the energy flux across the 
horizon inside the funnel gives a very similar number, 
Ebz — 0.42 x 10 51 ergs -1 , which is also in excellent 
agreement with eq.l. Thus, the explosion is undoubtedly 
powered by the Blandford-Znajek mechanism. Finally, 
this model seems to be marginal as any further notice- 
able increase of accretion rate or decrease of magnetic 
field strength prevents explosion. 

Although initially the BH jets are Poynting-dominated 
their electromagnetic energy is gradually converted into 
the energy of matter. To some extent this may reflect 
real physical processes (TT| [27l |28] |29) but numerical 
diffusion and resistivity are likely to play a significant 
part too. In any case, the total energy is conserved in 
the simulations and we do not expect these numerical 
problems to effect the dynamics of blast. 

In order to investigate further the potential of mag- 
netic breaking of accretion disk itself we carried out 
simulations with non-rotating BHs, thus eliminating the 
BZ-effect. The results show that the accretion disk is 
formed much later when almost all magnetic flux is al- 
ready accreted onto the BH and it begins to "swallow" 
the closed magnetic field lines. As in other solutions, 
we observe lift-off of strong accretion shock, develop- 
ment of disk corona, and amplification azimuthal mag- 
netic field. However, the explosion does not take place 
even under the most favourable conditions (the highest 
initial magnetic field and the lowest mass accretion rate) 
in our parameter space. This result provides an additional 
support to our conclusion on the crucial role of the BZ- 



mechanism. 

Our results seem to fit the following condition for 
magnetically driven explosion 

%i a5 *104,9 77- % > 0-1, (2) 

(1 + vl -^ ) 

where fl ,io = #o/10 10 G, r m , 9 r m /10 9 cm, M .i 
M/O.lM s _1 , which is obtained by comparing the max- 
imum power of the BZ-mechanism (see eq.l) and the 
power of accretion flow. Further investigation is needed 
to verify this condition for intermediate values of the spin 
parameter. For Mq.i > 0.5 the power of neutrino annihi- 
lation becomes substantial and this may lead to a much 
less demanding criterion. 

DISCUSSION AND CONCLUSIONS 

Our results suggest that magnetic fields may play a cru- 
cial role in powering GRB jets and associated hyper- 
novae not only in the magnetar but also in the collap- 
sar model. In the latter case the main source of energy is 
the rotation of BH and it is released via the Blandford- 
Znajek mechanism. The rate of energy release seen in 
the simulations, E = 0.2-1- 16 x 10 51 ergs _1 , is very high 
and can easily explain the energetics of hypernovae, E ~ 
10 52 erg. The fact that the rotational energy of maximally 
rotating BH is much higher, E rot ~ few x 10 53 erg, sug- 
gests a self-regulating process in which the rate of energy 
release significantly reduces after the first few seconds, 
presumably due to the lower mass supply into the ac- 
cretion disk when the stellar collapse is reversed by the 
blast, the escape of magnetic flux from the black hole, 
or its annihilation due to the accretion of magnetic field 
lines of opposite polarity. This idea is supported by the 
lower observed energetics of GRB jets. With the typical 
total energy of 10 51 erg and duration of 10 -r 100s their 
power is only 0.1 -=-0.01 x 10 51 ergs _1 EllED. Alterna- 
tively, the hypernovae blast wave could be caused mainly 
by the neutrino heating during the short initial period of 
high accretion rate. In this case much weaker magnetic 
field will be required to drive GRB jets into the cavity 
created by the blast. 
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